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Abstract 

The rise in linear entropy of a subsystem in the A^-atom Jaynes- 
Cummings model is shown to be strongly influenced by the shape of 
the classical orbits of the underlying classical phase space: we find 
a one-to-one correspondence between maxima (minima) of the linear 
entropy and maxima (minima) of the expectation value of atomic ex- 
citation Jz ■ Since the expectation value of this operator can be viewed 
as related to the orbit radius in the classical phase space projection 
associated to the atomic degree of freedom, the proximity of the quan- 
tum wave packet to this atomic phase space borderline produces a 
maximum rate of entanglement. The consequence of this fact for ini- 
tial conditions centered at periodic orbits in regular regions is a clear 
periodic recoherence. For chaotic situations the same phenomenon 
(proximity of the atomic phase space borderline) is in general respon- 
sible for oscillations in the entanglement properties. 
PACS numbers: 05.45.Mt, 32.80.Qk 
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I INTRODUCTION 



The importance of studying in detail the decoherence process is twofold. 
Firstly it may be viewed as a key to the understanding of some of the strik- 
ing differences between the quantum and classical description of the world 
such as "the non existence at the classical level of the majority of states al- 
lowed by quantum mechanics" The decoherence process is believed to be 
the agent which eliminates interference between two or more macroscopically 
separated localized states [^. Secondly, given the impressive technological 
advances in several experimental areas (Quantum Optics, Condensed Mat- 
ter and Atomic Physics, etc.), it is nowadays possible to realize a system 
of two interacting degrees of freedom and watch the time evolution of the 
corresponding entanglement process Q]. It is therefore also of importance to 
understand the entanglement process in simple Hamiltonian systems. Two 
degrees of freedom Hamiltonian systems often present a very rich dynamics 
which in many cases is not yet completely understood from a general point 
of view. In particular, if the interaction is nonlinear the system may present 
chaotic behavior in the classical limit. The consequences of this fact to the 
quantum dynamics is yet an unsettled issue. A step in this direction was 
taken a few years ago, as it was conjectured that "the rate of entropy pro- 
duction can be used as an intrinsically quantum test of the chaotic versus 
regular nature of the evolution" [Q. The idea has been tested in some mod- 
els 1^, ^. More specifically, in the context of the A^-atom Jaynes-Cummings 
model, the reduced density linear entropy (or idempotency defect) has been 
used as a measure of the entanglement of the quantum subsystems. For a 
given classical energy, initial conditions for the quantum states are prepared 
as coherent wave packets centered at regular and chaotic regions of the clas- 
sical phase space. For short times, a fast increase in decoherence for chaotic 
initial conditions is found when compared to regular ones. Typically the 
linear entropy in this model rises from zero to a plateau. 

In the present contribution we show that this rise is strongly influenced 
by the shape of the atomic projection of classical orbits. There is a clear 
correlation between the rate of increase of linear entropy and the increase 
of the expectation value of the atomic excitation {Jz){t) for short enough 
times. In classical terms, this can be visualized as follows: {Jz){t) is a mea- 
sure of the instantaneous radius of the projection of the trajectory in the 
atomic phase space. This atomic phase space is limited, the maximum ra- 
dius corresponds to v^4J = v^2iV, meaning therefore that the maximum rate 
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of entanglement can be directly associated to the proximity of this border- 
line. If the quantum evolution is such that the initial wave packet is centered 
at a classically regular region of the atomic phase space, in particular on a 
(non circular!) periodic orbit, the phenomenon is most conspicuous: as a 
function of time, the wave packet exhibits partial recoherences, well marked 
oscillations in the linear entropy superimposed to a steady increase, which 
can be immediately associated to smallest and largest distances from the 
classical border (largest or smallest {Jz){t)). If, on the other hand, the initial 
wave packet is centered on a classically chaotic region of the phase space, 
the recoherence phenomenon is of course not necessarily periodic in general. 
However we still find a connection between maxima and minima of the rate 
of linear entropy and ( J^) (t). As to be expected, this simple relation becomes 
less clear as quantum effects set in; in particular, for the chaotic region this 
happens much faster than for the regular one in a mixed phase space. Such 
effects are most marked at sufficiently long times, when a plateau in the linear 
entropy is reached. We show that at these times the spin projection of the 
wave packet becomes totally delocalized in phase space. This takes place for 
initial conditions centered at regular and chaotic regions, the regular cases 
taking longer times to achieve this delocalization. 

In section II we present the model and the technical tools to be used in 
the following sections to analyze the entanglement property, connecting it to 
the classical regularity and chaos in a mixed regimen. Section III is devoted 
to the presentation of analysis and results in regular regions, whereas section 
IV deals with chaotic regions. A summary is given in section V. 



II THE MODEL 

We consider the iV-atom Jaynes-Cummings model whose Hamiltonian is 
given by 

G G' 

H = huoQa^a + eJz + "77^ {aJj^ + a+ + {a^Jj^ + aJ_^ (1) 



where the first term corresponds to the energy of the free single-mode quan- 
tized field with frequency ujq, described by the creation (annihilation) oper- 
ators a+ (a); the second term corresponds to the energy of the N = 2J two- 
level atoms with energy separation he, and the operators Jz, J± = Jx ± iJy 
are the usual angular momentum ones corresponding to the group SU{2). 
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These operators satisfy the well known commutation relations: 



a,a~^ = 1, 
[Jj, Jj] = ihJk 

with the indices k) forming any cyclic permutation of {x,y, z). The last 
two terms in eq. ([1|) represent the interaction energy between the atomic 
system and the single-mode field. 

The above Hamiltonian is a generalized version of the usual iV-atom 
Jaynes-Cummings model to which we have added a Counter-Rotating Wave 
(CRW) term, with coupling constant G . The so-called Rotating Wave Ap- 
proximation (RWA) has been much in use for the proposals of mesoscopic 
superpositions of collective atomic states R P[ and in the generation of mul- 



tiparticle entangled states in cavity QED experiments |]T0 . 

The iV-atom JCM can also be realized in systems where trapped ions 
interact with laser fields. In that situation the bosonic operators a and a+ 
describe the vibrational motion of the ions in the trap. Recently, reversible 
entanglement of ions has been proposed and experimental realizations 



have been reported in the cases of two and four ions |jT2| , with possible exten- 



sions to several ions |T3|. In the trapped- ion system there is the possibility 
of generating both RW and CRW types of interaction by means of either 
a bichromatic laser excitation — with both G and G different from zero 
— or appropriate excitation modes fl^. In the latter case, there would be 
the possibihty to test the effects of RW (G' = 0) and CRW {G = 0) terms 
separately or together. 

In our model we also add the supposition that the two coupling constants 
can be independently varied. When one of the coupling constants G or G 
is set to zero, we have an integrable system. For G = [G = 0) there is 
an additional conserved quantity, namely the total excitation P = J^/h + 
a^a (the relative excitation P = Jz/h — a~^a). The system is otherwise 
nonintegrable and known to exhibit chaotic behavior in the classical limit 



[jT5[. In both cases, the quantum dynamics will produce entanglement, due 
to the coupling terms. This essentially quantum property can be quantified, 
e.g. by means of the linear entropy (or idempotency defect) associated to the 
reduced density matrices of the atomic or field subsystems (from this point 
on we shall set h = 1), 

5(t) = l-Tr,(p,(t)2), (2) 
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where 



p,{t) = Tr, mt)) m)\) (3) 

(i and j stand for the atomic and field subsystems, i 7^ j). The quantities 
expressed in eqs. (||) and (H) are to be calculated using quantum states 
evolving in time under the action of Hamiltonian (|l]). To this end, the initial 
conditions chosen in the present study are such that 

\^{0)) = \w) (g) \v) = \wv) (4) 

where \w) ( \v) ) are atomic (field) coherent states given by |I6 



\w) = (l+ww)-^ e'"-^+\J,-J) (5) 
It;) = e-™/2e^''^ |0) (6) 

with 

Pa + iQa 

W = I (7) 

/4 J - {Pl + ql) 



^ = ^ (p/ + nf) ; (8) 

|J, — J) being the state with spin J and Jz = —J, |0) being the harmonic 
oscillator ground state, and Pa,(la,Pf,(lf describing the phase space of the 
system into consideration. Generation of atomic coherent states in A^-atom 
systems has been proposed in and analyzed in 0. 

The classical Hamiltonian corresponding to eq. (|T]) can be obtained by a 
standard procedure as 

H{v,v* ,w,w*) = {wv \H\wv) . (9) 
7i(f , V*, w, w*) can be rewritten in terms of the phase space variables, reading 

T-f-iQa, Pa, Qf, Pf) = ^ (Pj + (if) + ^ (Pa + ?f - 2J) + 



1 - (G+PaPf + G^q^qf) , (10) 



where G+ = G±G'. 
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The corresponding equations of motion are then given by 



^« ~ ~~a — ~ -^Pa - ^+Pf\ — h —?= / 

(11) 

It is important to note the restriction in energy to be shared in the clas- 
sical counterpart of the atomic degree of freedom, i.e., 

pj + qJ < 4J. (12) 

This represents a fundamental difference between these two degrees of free- 
dom. From the quantum point of view, the atomic degree of freedom is 
associated to a finite dimensional Hilbert space whereas the harmonic oscil- 
lator to a Hilbert space of infinite dimension. 

Except when otherwise noted, in all cases to be presented here the pa- 
rameter values we chose are J = 9/2, G = 0.5 and G = 0.2, and e = ujq = 1 



(as shown in Ref. [0, for these parameter values the nearest neighbor dis- 
tribution of the energy level spacings is of a GOE type). For the classical 
Hamiltonian in eq. (|10]) we chose an energy value E = 8.5 for which the 
Poincare section (gj = 0) is shown in Fig. la. Note the presence of two is- 
lands, one of them of considerable relative size. In the center of such islands 
are located the two shortest stable periodic orbits whose atomic projections 
are shown in Fig. lb. The symbols in the Poincare section (Fig. la) repre- 
sent the center of the initial quantum wave packets whose time evolution we 
will show in the next section. For the sake of comparison we also study the 
classical energy = 35 whose Poincare section is shown in Fig. 2a. The size 
of the many regular islands is comparatively much smaller. Both chosen en- 
ergies are larger than the limit given by eq. ([T^), therefore the entire atomic 
phase space is energetically available for all initial conditions in the Poincare 
sections. Moreover, the fact that the energies are large enough together with 



the limitation imposed by eq. (12) have dramatic consequences on the form 
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of the projections of the periodic orbits onto the atomic phase space. One 
can see an example of this in Fig. 2b (this particular point has been explored 
inRefs. |Ti). 



Ill ENTANGLEMENT DYNAMICS IN 
REGULAR REGIONS 

The linear entropy 6{t) corresponding to the reduced density matrix for the 
atomic subsystem is shown in Fig. 3a for initial conditions of the type given 
in eq. where the center of these wave packets are on the periodic orbits 
indicated in Fig. la. Note that the linear entropy not only increases but 
exhibits a well marked oscillatory behavior, indicating that the system re- 
covers coherence in a periodic way, for sufficiently short times. Interestingly 
enough, however, the period of recoherences observed in the linear entropy 
do not correspond to the period r of the periodic orbit in question. It is 
smaller, roughly r/2. This indicates that another property of these orbits 
may be playing a role. This property is the shape of the spin (atomic) 
projection of the orbit. Note the strong correlation between the maximum 
growth in linear entropy and the closest approach to the classical atomic 
phase space border. This happens in the cases shown here and in all cases 
we analyzed, for short enough times. A natural question arises as to the 
representation dependence of our explanation. In order to clarify this point 
we use a representation independent quantity which is intimately connected 
to the radius of the orbit projection, the expectation value of the operator 
J2 pO|]. In Fig. 3a we also show ^ and Tr{JzPa{t))- Note that the first four 



maxima in {Jz){t) correspond to a very good accuracy to maxima in ^, indi- 
cating that whenever the "radius" is maximum, a maximum growth in linear 
entropy is found. It is in this sense that we conclude that the recoherences, 
purity gains found in the time evolution of quantum coherent wave packets 
initially centered at the regular region, are related to the proximity of the 
classical atomic phase space border. Of course when the initial wave packet 
is centered on a periodic orbit located in a smaller island, as is the case of the 
second orbit in Fig. lb, the chaotic vicinities also play an important role and 
the established connection between ^ and Tr{JzPa{t)) is less pronounced, as 
shown in Fig. 3b. Due to the shape of the orbit, in this case the proximity 
to the center and hence the minima of both functions are more clearly seen 
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to be connected. We next consider two more cases which independently cor- 
roborate these findings. Consider first the Poincare section in Fig. 2a and 
the initial condition for the quantum evolution centered at the circle marked 
in the section, and the spin projection of the periodic orbit shown in Fig. 
2b. In this case we note that the amount of time spent in the vicinities of 
the phase space border (performing the "loop" in Fig. 2b) is large compared 
to the period of the orbit. Note that during this time, ^ is positive and the 
linear entropy grows (see Fig. 3c). Also, for the subsequent times where the 
orbit quickly crosses from one "loop" to the other, it approaches the center 
of the atomic phase space and correspondingly ^ becomes negative and the 
linear entropy decreases. Since the time of growth is much larger than the 
one of decreasing, the linear entropy reaches the plateau faster. 

On the other hand, in the integrable case G = 0, the spin (atomic) 
projection of the periodic orbits are circular. In this case, if the argument is 
valid, one should expect no oscillations in 6{t). This is in fact the case, as 
shown in Fig. 4 where one sees, moreover, that it takes a longer time for the 
linear entropy to reach the plateau, due to the considerable (and constant) 
distance between the orbit and the atomic phase space border. 

IV ENTANGLEMENT DYNAMICS IN 
THE CHAOTIC REGION 

As discussed in Ref. ^ the entanglement in the chaotic region is in gen- 
eral faster than in the regular region. In Fig. 5a we show the idempotency 
defect (linear entropy) for the initial conditions of the form (^) centered at 
the triangles marked on the Poincare section of Fig. la, and in Fig. 5b the 
projection of the orbits for the specific cases of the initial conditions labeled 
Ci and C2 (see figure captions). In these cases the linear entropy does not ex- 
hibit periodic oscillations, albeit oscillations are still present. Although they 
are not directly related to a periodic orbit, one can explain the oscillations 
very much along the lines followed in the regular case. If we take a classical 
trajectory associated with the center of the wave packet and follow its pro- 
jection on the atomic phase space, we notice again that the proximity of this 
phase space border plays a decisive role on the decoherence process. This is 
illustrated in Figs. 5c and 5d for the initial conditions ci and C2- For the case 
C2, decoherence is retarded due to the fact that the wave packet moves ini- 
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tially towards the center of the phase space, following the classical trajectory 
shown in Fig. 5b. Note that as soon as the classical trajectory approaches 
the border there is a corresponding increase in the entanglement rate ^ in 
Fig. 5d. Thus the linear entropy increases. In terms of the representation 
independent quantity {Jz){t), the same features prevail. The remaining ini- 
tial condition C3 for which 6{t) has been evaluated in Fig. 5a behaves in 
qualitatively analogous way, as for instance the trajectory labeled Ci shown 
in Fig. 5b, and will not be shown here. As to be expected, in the chaotic 
region the relation between the motion of the quantum wave packet and the 
classical trajectory associated to its center is rapidly overwhelmed by other 
effects such as the rapid spreading of the quantum wave function. In fact, in 
Fig. 6 we show the time evolution of the spin Husimi distribution associated 
with the initial condition C2 in Fig. 5. This figure illustrates the given argu- 
ment. Moreover, we see that for times when the plateau in 6{t) is reached, 
the spin Husimi distribution is totally delocalized in the corresponding phase 
space. Finally we note that this is not a characteristic of the chaotic initial 
condition; the wave packets centered at regular initial conditions also become 
spread all over the spin phase space by the time the linear entropy reaches 
the plateau. 



V SUMMARY 

The present work has been devoted to a detailed analysis of the dynam- 
ics of the process of entanglement in the chaotic A^-atom Jaynes-Cummings 
model. Previous work on the model ^ pointed out general features of the 
process, mainly focusing on differences between chaotic and regular regimes 
(unstable and stable). The entanglement process of the model is however 
very rich and much more can be learned from specific features of the linear 
entropy such as its oscillations (recoherences in time, periodic or not). We 
have found an intimate connection between these recoherences and classical 
orbits. Specifically, we have pointed out the very special role played by the 
morphology of the spin-projected classical orbits. This is in accordance with 
the known fact that the decoherence properties are dictated by the smaller 



subsystem [22 



Considering the present state of experimental proposals and realizations 
of atomic coherent states in the A^-atom JCM, it seems that the integrable 
RW case {G = 0) is more likely to be realizable in a practical implementation 
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of the model considered in this work. Although not shown here, the RW case 
was checked for regular initial conditions on tori in phase space and showed 
similar recoherence phenomena, as described in the previous sections. 

Despite the fact that these findings are model dependent, we believe them 
to be typical of systems involving two degrees of freedom, one of them having 
a Hilbert space with dimension much greater than the dimension of the other, 
and whose classical hmit is chaotic (soft chaos). 
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Figure Captions 

Figure 1: (a) Poincare section given by Qf — 0.0 and pf > 0.0 for 
the spin degree of freedom with J = 9/2 in the nonintegrable 
{G = 0.5 and G' = 0.2) and resonant {e = wq = 1) case — 
these values wiU be the same in all subsequent figures except 
when otherwise noted. Here we show the section for the energy 
E — 8.5. The marks represent the various choices for the center of 
the coherent states: circles for regular initial conditions (IC) and 
triangles for chaotic ones, (b) Spin projection of the two shortest 
stable periodic orbits indicated by circles in (a) (the border line 
is also shown). Initial conditions at the surface of section are 
as follows: (g„ = 0.0,p„ = 2.261, g/ = 0.0,p/ = 3.423276) for 
the first (elliptic) orbit with period r = 4.89; (g^ = 0.0, = 
—3.577, Qf — 0.0, p/ = 5.221656) for the second orbit with period 
r = 7.45. 

Figure 2: (a)Poincare section for the energy E — 35. The circle repre- 
sents the center of the coherent state for one regular initial condi- 
tion, (b) Spin projection of the shortest stable periodic orbit, in- 
dicated by the circle in (a). Initial condition at the surface of sec- 
tion is as follows: (g„ = 0.0, p„ = 1.4175, qf = 0.0, p/ = 7.888904) 
with period r = 5.82. 

Figure 3: Linear entropy 6{t) (continuous lines), its derivative 

(dashed lines), and the normalized expectation value {Jz){t)/J 
(dot-dashed lines) as functions of time, (a) for the first regular 
initial condition of Fig. lb. Note the correspondence between the 
maxima and minima of and {Jz){t)/J. (b) for the second 
regular initial condition of Fig. lb. In this case, only the minima 
of and {Jz){t)/ J have a clear correspondence due to the pe- 
culiar form of the orbit (see second figure of Fig. lb), (c) Regular 
initial condition as marked in Fig. 2a and its spin-projected orbit 
shown in Fig. 2b; during the interval of time corresponding to 
the rise of 5{t) there is an approximate correspondence between 
the maxima and minima of and those of {Jz){t)/ J. 
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Figure 4: For the energy E — 8.5 in the integrable (G = 0.5 and G' — 
0.0) and resonant {e = Wq = 1) case, the plot of hnear entropy 
6{t) for initial condition on a periodic orbit which is circular (see 
inset) {qa = 0.0,p„ = 2.47675, g/ = 0.0, pf = 3.563642). Notice 
that the curve shows practically no oscillations as compared to 
Fig. 3a. 

Figure 5: (a) Linear entropy 6{t) for the chaotic initial condi- 
tions corresponding to the triangles shown in Fig. la: Ci 
{qa = -4.0,p„ = 0.0, qf = 0.0,p/ = 3.162278) for the continu- 
ous line; C2 {qa = 1.57, = — 2.0,g/ = 0.0, pj = 5.680465) for the 
dashed line; and C3 (g„ = 3.0, = 2.0, g/ = 0.0, pf = 2.942413) 
for the dot-dashed line, (b) Spin projections of the orbits ci and 
C2. Various times are indicated along the orbits showing that the 
initial conditions are such that the trajectory C2 in spin phase 
space is launched inward as opposite to the case Ci for instance, 
(c) S{t) (continuous line), ^^^^ (dashed line), and {Jz){t)/J (dot- 
dashed line) for ci. (d) 6{t) (continuous line), (dashed line), 
and (Jz) {t) I J (dot-dashed hue) for c^. During the interval of time 
corresponding to the rise of 5(i) there is an approximate corre- 
spondence between the maxima and minima of and {Jz){t) / J 
but for shorter time than it happens for the regular cases shown 
in Figs. 3a and 3b. 

Figure 6: Spin Husimi distribution of the quantum coherent state ini- 
tially centered at condition C2 in Fig. 5. Snapshots are taken at 
time values t = (upper left), t = 1 (upper right), t = 4 (down 
left), and t = 25. 
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